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The method developed in [L. P. Gor'kov and L. P. Pitaevskii, Sov. Phys. Dokl. 8, 788 (1964); C. 
Herring and M. Flicker, Phys. Rev. 134, A362 (1964)] to calculate the asymptotic form of exchange 
interactions between hydrogen atoms in the ground state is extended to excited states. The approach 
is then applied to shallow centers in semiconductors. The problem of the asymptotic dependence of 
the exchange interactions in semiconductors is complicated by the multiple degeneracy of the ground 
state of an impurity (donor or acceptor) center in valley or band indices, crystalline anisotropy and 
strong spin-orbital interactions, especially for acceptor centers in III-V and II- VI groups semicon- 
ductors. Properties of two coupled centers in the dilute limit can be accessed experimentally, and 
the knowledge of the exact asymptotic expressions, in addition to being of fundamental interest, 
must be very helpful for numerical calculations and for interpolation of exchange forces in the case of 
intermediate concentrations. Our main conclusion concerns the sign of the magnetic interaction — 
the ground state of a pair is always non-magnetic. Behavior of the exchange interactions in applied 
magnetic fields is also discussed. 
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I. INTRODUCTION 

Ferromagnetism in diluted III-V semiconductors has 
recently attracted considerable attention for technologi- 
cally appealing possibility of integration of spin degrees 
of freedom into semiconductor devices. In the most inten- 
sively studied compound of this group Gai-^Mn^As the 
ferromagnetic transition temperatures in thin annealed 
epitaxially grown layers as high as Tc > 150K for x ~ 8% 
have been reported^. 

The transition metal impurities in these materials sub- 
stitute for the atoms of the III group and act as both 
localized spins (S = 5/2 for Mn due to half-filled d-shell) 
and as acceptors doping the system with one hole. How- 
ever, the system is highly compensated presumably owing 
to double donor antisites (V group elements substituted 
for III group elements) and metastable interstitial tran- 
sition metal impurities. 

A major not fully resolved issue is the mechanism gov- 
erning ferromagnetism in these compounds. Despite high 
disorder in these non-equilibrium epitaxially grown films, 
RKKY interaction seems to be the most plausible cause 
of Mn spin-spin interaction^. 

Meanwhile, we want to indicate that magnetic in- 
teractions in semiconductors are not well understood 
yet even in a much simpler case of electrons localized 
on non-magnetic impurities. Whereas two hydrogen 
atoms at large distances are known to interact anti- 
ferromagnetically, the sign of the interaction between 
donors/acceptors is not obvious because of the band 
degeneracies as well as strong spin-orbital interaction 
present. 

To the best of our knowledge, exchange interaction be- 
tween donors was considered only in Refi^ in the Heitler- 
London approximation which is known to be incorrect 
at large separations^^. Besides, in Refi^ an isotropic 



hydrogen-like model was adopted in contrast to real semi- 
conductors that have strong anisotropy of conduction 
band valleys. In Section ITTT1 below we compare our re- 
sults for donors to those of RefA 

The purpose of the present paper is to obtain asymp- 
totically exact formulas for the exchange interaction be- 
tween two shallow centers at large separations. Although 
our results are pertinent to low impurity concentrations, 
this information would be necessary for interpolating re- 
sults of numerical calculations to shorter inter-impurity 
separations — higher impurity concentrations. 

The Heitler-London scheme, which considers exchange 
with the alien atom as perturbation, is inapplicable^ to 
electron exchange at large separations R, because the 
tunneling takes place through a region where interac- 
tions with both atomic residues and between the elec- 
trons themselves are comparable. As a matter of fact, the 
Heitler-London expression for the exchange between two 
hydrogen atoms incorrectly changes sign at large inter- 
atom separations. The correct approach to the asymp- 
totic exchange, that accurately accounted for the interac- 
tion between the electrons and ionic residues during the 
exchange process, was proposed in***&. The magnitude of 
the exchange was expressed as an integral of the proba- 
bility current over the 5-dimensional median hyperplane 
separating two regions of electron localization in their 
6-dimensional coordinate space. 

Before proceeding to the corresponding problem in 
semiconductors we list a number of extensions of the 
method*^ to other problems in atomic physics. In£ a gen- 
eralization was made for two atoms with non-equal bind- 
ing energies of electrons in s-states. It turned out, how- 
ever, that in order for the method to be applicable, the 
two energies should be close, which severely restricted the 
applicable range of R (see the discussion in&). In Ref. 9 
exchange integrals between degenerate s- and p- states 
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have been calculated starting from Wigner-Witmer clas- 
sification of the molecular states (see alsoAfi) . Exchange 
integrals between higher orbital momentum I states were 
calculated ir*ii. In papers^ the method*^ was applied to 
get interpolated formulas down to small R and to study 
two-dimensional electrons localized on out-of-plane im- 
purities. But to the best of our knowledge, multiple 
degeneracy of donor or acceptor levels in real semicon- 
ductors has never been addressed. 

In this paper we extend the median-plane method to 
asymptotic exchange between shallow impurities in semi- 
conductors with localized carriers in degenerate and non- 
spherically symmetric bands. We also study how spin- 
orbital interaction and external magnetic field influence 
the exchange. 

As is well known, the basis of description of shal- 
low impurity centers in semiconductors is the effective 
mass metho d 1 ^' 14 . Shallow ceneters have small ioniza- 
tion energy compared to the forbidden gap, and, con- 
sequently, the localization scale much larger than the 
lattice constant. Hence for shallow centers the poten- 
tial binding carriers to the oppositely charged impurity 
ion may be deemed of approximately Coulomb form: 
V(r) — —e 2 /nr, where k is the permittivity. Following 
standard Kohn-Luttinger formalis m 1 ^ 14 , the wave func- 
tion may then be expanded in many-electron Bloch func- 
tions ip( a > (r) at the extremum of the free band 

^ a \r) = <^ Q )(r)V> (Q) (r)> (1) 

where a = 1, . . . , F, and T is the band degeneracy at the 
extremum, and the envelope function ip a (r) satisfies a 
Schrodinger equation 

E<p a = -\DgdidjtpP + V{v)<p a (2) 

with the effective mass tensor determined by the 

crystal symmetry. 

We have generalized the approach^ to the form of 
a secular equation, and in Section [H] of the paper we 
demonstrate the method on the example of two atoms 
with degenerate states. The atomic problem has much 
in common with the case of two acceptors. In Section 
IIIII we treat shallow donors in III-V (and also group IV) 
semiconductors where the conduction band has several 
equivalent minima. Section IWl addresses the question of 
asymmetric exchange in III-V semiconductors. Section 
Ivl deals with acceptor impurities localizing holes from 
degenerate valence bands. In Section Ivfl the dependence 
of the exchange on the applied magnetic field is studied. 
In the closing section we summarize our results. 

II. TWO HYDROGEN-LIKE ATOMS 
SECULAR EQUATION APPROACH 

Consider two identical atoms A and B located at 
— R/2 = —RC/2 and R/2 respectively, each having one 
valent electron outside of closed shells. Let an isolated 



atom have degenerate states tp v ( r ) with binding energy 
—a 2 /2. Here v are the quantum numbers describing the 
state. In central-field approximation v stands for a set 
of principal quantum number n, orbital momentum I, its 
projection m and spin projection s. Henceforth we use 
atomic units, i.e. measure lengths in units of the Bohr ra- 
dius ae = h 2 /m e e 2 , energies in units of 2Ry = m e e 4 /H 2 . 
At large distances 

¥v(r) « C u Yrmri- l e- ar , (3) 

where Y"j m (f ) is a spherical function, and C v is the asymp- 
totic coefficient. 

The two-electron wave function ^"^(ri,^) in the ab- 
sence of spin-orbital interaction can be factorized into 
the product of spin x and coordinate ip(r 1 ,r 2 ) parts, 
where the former may be either a singlet S — or a 
triplet 5=1, and the latter is an eigenfunction of the 
Hamiltonian 

H = -\d\-\d\-\lrxK-\lr-m (4) 
- 1/raA - 1/riB + l/r u + 1/R. 

Here S = Si + S2 is the total spin and Tia = t% + R/2, 
Tib =ri— R/2, etc. 

In <Q} the first row is the Hamiltonian of two isolated 
atoms, while the terms in the second row result from 
bringing the two atoms to a finite separation. These 
terms are taken as perturbation in the Heitler-London 
approximation. The correct treatment of these terms has 
to be accomplished in two stages. 

First, if the possibility for the electrons to tunnel be- 
tween the two potential wells is neglected, this part of 
the Hamiltonian may be expanded in multipole moments 
over inverse powers of R. Perturbation theory calcu- 
lations then yield van der Waals corrections to the en- 
ergy of each electron. We denote the energy of two elec- 
trons with van der Waals interaction taken into account 
as Eq. Corresponding corrections to the direct product 
<Pv A ( r i)¥VB ( r 2) of isolated-atom wave functions retain 
electron localization. We denote these perturbed wave 
functions when the first electron is localized on ion A in 
the state v\ and the second — on ion B in the state 

aS IpL/A I^B • 

Second, overlap of the wave functions of electrons © 
on the two atoms enables them to exchange position 
by tunneling. This lifts the original degeneracy in en- 
ergies of the isolated-atom states. In contrast to van 
der Waals interactions, the resulting exchange splitting, 
which exponentially decreases with distance, cannot be 
calculated perturbatively. One can use symmetry reason- 
ing to construct the correct two-electron wave functions 
from the localized wave functions ip VAVB and interchanged 
P^vavb = VVb^a- Then the energy splitting is found by 
substituting those combinations into Hamiltonian (@J). 

In the simplest case when atomic states are each only 
doubly degenerate in spin projection s, the correct two- 
electron functions are symmetric and antisymmetric com- 
binations of the localized wave functions^. As total 



3 



fermionic wave function must be antisymmetric in parti- 
cle permutation ip a(3 (ri, r 2 ) = — ^ a (r 2 , i"i), the singlet 
state corresponds to a symmetric coordinate wave func- 
tion and the triplet state — to an antisymmetric one. 
Hence the energy splitting between triplet and singlet 
states is given by 

E ex = -2J{R)( Sl s 2 + l/4), (5) 

where J(R) oc e ~ 2aR is found by substituting symmet- 
ric and antisymmetric combinations into Hamiltonian <@} 
(see below). 

However, if degeneracy of atomic states is higher, the 
choice of correct combination of two-electron wave func- 
tions becomes more complicated. Higher orbital momen- 
tum I states of two atoms are degenerate in its projec- 
tion m. In this case, the Wigner-Witmer classification of 
molecular states in the absolute value of the total an- 
gular momentum A (and also parity for homo-atomic 
molecules) can be used to identify the correct combina- 
tions of localized wave functions^. The combination sat- 
isfying the required symmetry constrains, however, can 
sometimes be not unique. Or else, such a classification in 
cases other than atomic may be complicated or absent. 

We argue that one can simplify the procedure and find 
the correct coefficients in the expansion 

1p= ( C ^a.b^a./b+c1aIb^a.b) (6) 

as well as the energy level splitting in an approach simi- 
lar to the secular equation of perturbation theory of de- 
generate level. To this end one substitutes © into two- 
electron wave equation with the Hamiltonian J3J and uses 

In the next few paragraphs we derive an analogue for 
secular equation in this approach. Let us denote the set 
of indices va^b by a single letter n. We start from the 
matrix continuity equation ihdt4>iipk = —ifrdjiki which 
is a direct consequence of the Schrodinger equation and 
where the probability current density 

3ik = ~W>ifyk-4>ktol>i)- (7) 
2m 

Writing it for ipi — ip n and ipk = ip from @ and inte- 
grating over the region £i < £2, we obtain 

(E - E Q ) f C^n,2 - -iH f j„^dS. (8) 

Here dS = (£ l7 — C 2 )^C^ 4 Pi, 2 i where p 12 are the trans- 
verse radii-vectors of the electrons, is the element of 
the median 5-dimensional hyperplane £ = £1 = C2 that 
bounds the integration region. 

Functions ip n describe electrons localized on different 
ions, therefore the integral in the left hand side of Eq. 
{HJ equals c„ up to an exponentially small quantity which 



we neglect as being of the same smallness as the energy 
splitting E — Eq itself. From JJJ we then obtain 

(E - E a )c n =^2 / [(^m% 2 C - K 5 Ci2^m)c OT 

m "'Cl=C2 

+ (P^ m d Cl2 r n -r n dc 12 P^ m )c^]dCd 4 p la , (9) 

where £ = ±(£1 + £ 2 ) and £ X2 = £1 - £ 2 . 

Differentiating the functions ip m and Pi/ 1 ™ in this ex- 
pression one may retain only derivatives from the rapidly 
decaying exponents ip m ~ e - Q (-R+Ci2) ; so d^ 12 ip m « 

—aijj m , and d^ 12 Pip m s» aPip m . In this approximation 
the expression in the first parentheses vanishes and we 
arrive at 

{E-Eo)cn = J2 J ™ c m } ' ( 10 ) 

m 

with the exchange integral matrix elements given by 

J nm = -2a / tl)* n Pi) m dC,d 4 p 1 2 . (11) 

" / Cl=C2 

Likewise, starting from the permuted states Pip n we 
integrate over the region £1 > £2- Combining the result 
with ijlOjl . we get 

(*LV )(&)-<* -*>(£) <») 

(the Einstein convention on summation is implied). Since 

coordinate wave functions of the triplet states are odd in 

(p) 

electron interchange (c„ = — c„ ) and those of the singlet 
states are even (c„ — c n ), eigenvalue equation simplifies 
to 

JnrnC m = ±(£7 — E )c n , (13) 

where sign plus corresponds to singlet (even) solutions, 
sign minus to triplet (odd) ones. We conclude that the 
splitting of originally degenerate state \va) ® I^b) due to 
exchange must be found as eigenvalues of the exchange 
integral matrix J VAUB y A v' and the coefficients in combi- 
nations 10 of localized wave functions as corresponding 
eigenvectors. 

The rest of the Section is devoted to the details of 
evaluation of exchange integrals in atomic case, which 
will come useful for treatment of acceptors below. To 
calculate J UAVB y v ^ one needs to know the localized wave 
functions ip VAUB . We are only interested in asymptotic 
dependence of exchange at large R, hence we will find 
VVa^b m principal order in R. As m&&, we seek i(j vaub in 
the form 

^ A u B (ri,r 2 ) = x^i^2)Vv a {yia)^u^2b), (14) 

where x is a smooth function varying on the scale of 
R. Substituting 11411 into the Schrodinger equation 
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Hip VA v B = Eoip VAVB , where E = -a 2 + 0(i?~ 6 ) and 
keeping only terms of order R^ 1 (and thus neglecting 
the second derivatives of x) > we obtain 



(15) 

Introducing the new variables C = |(Ci + C2) and C12 = 
Ci - C2, P12 = I Pi - P 2 |, in which % - % = 2% 2 , the 
solution of this equation is 



x = c(Oi 2 ) 



+ P? 2 " Cl2 



(i?-Ci2-2C)(i?-Ci2+2C) 



-C12/2-R 



(16) 

where the integration constant C(C,pi2) must be deter- 
mined from the boundary condition % — > 1 when either 
n -> -R/2 or r 2 -> R/2. 

On the median plane Ci = C2 = C one thus obtains 



R 



P12 



R/2 + ICI V A/2 - ICI 



CI-H/2 

e h 



(17) 



We now substitute ifPIj) into the integrals in exchange 
elements (11 II . Expanding the asymptotics JSJl: e~ ar « 
e -"(ICI+p /2|CI) ) the integrals over transverse coordinates 
Pi 2 acquire the Gaussian form and converge rapidly on 
scale p ~ We thus need the asymptotics (pj only 

in the narrow cylinder p ~ V^R <C C ~ R, i- e - at 
small polar angles 9 from the internuclear axis C chosen 
as the quantization axis for orbital momentum. Since 
Y l m (e,ip)^pY i e^e imv /V^ when 6 < 1, where 



Pi 



(-)- 



2M|? 



'(21 + 1) 



q + H)! 

(I-H)!' 



(18) 



Eq. (pj may be rewritten as (9 » p/|Cl) 



^( r ) ~ _^pl'»l|ds- 1 -l m l e -«(ICI+P 2 /2|CI) e ^ ) ( 19 ) 



47T 

where A„ = CpJ" for atom A, and A„ = (-l) ! - m C^ m 
for atom B. The sign multiplier for atom B originates 
from the change of 9 to tt — 9. 

If all to = the exchange matrix elements J may be 
evaluated with the help of the substitution 

P=\{Pl+ P2)i Pl2=Pl-P2i (20) 

which yields! 



Jv 



'VAV B V A V B V A i-*B "A "B 

where 

™i 



A* i4J B A^^/(a)i2 r/2a - 1 e- 2ttfl ) (21) 



(22) 



For non-degenerate s-states p[J = l(see ill8|) ). Ao = Co 
and the energies of singlet and triplet states are 



E. tt ^EoTlCo^H^R 7 ^-^- 



2aR 



(23) 



The asymptotic parameter Co of the atomic function 
(GJ depends on the generally unknown form of the atomic 
potential V(r) at small r. For a hydrogen atom in an s- 
state Co = 2, a = 1, 1(1) « 0.511 and we regain that 
the ground symmetric singlet and excited antisymmetric 
triplet terms are separated bj£ 2 • 0.818i? 5/2 e~ 2fl . 

To find the hierarchy of exchange splittings in the gen- 
eral case, one needs expressions for J for other values of 
to given in the Appendix. 



III. EXCHANGE INTERACTION BETWEEN 
SHALLOW DONORS 

In Ge, Si and some zinc blende semiconductors (GaP, 
AlSb) several equivalent conduction band minima are lo- 
cated in equivalent points of the Brillouin zone, turning 
one into each other under cubic symmetry transforma- 
tions. In Ge these are four L-points of intersection of 
third order A-axes [111] with the Brillouin zone bound- 
ary. In Si and III-V semiconductors minima lie on the 
three fourth order A-axes [100] at the A-points (III-V) 
on the Brillouin zone boundary or close to them (Si). 
Electrons in the coaxial pair of valleys in Si are indistin- 
guishable, and such a pair may be considered as one with 
doubled density of states. 

Isoenergy surfaces in the vicinity of a minimum num- 
bered by the index u are ellipsoids of revolution around 
the symmetry axis d u , and the effective-mass tensor has 
uniaxial structure 



D„ 



dui duj 



) + m H 



(24) 



where are the transverse and longitudinal effective 

masses of electrons and d U i are the direction cosines of 
the axis d„ in a certain coordinate system. In atomic 
units defined now as ae = n/m±e and 2Ry = m±e /k , 
i.e. corresponding to the transverse mass: 



E*uij K^i; 



(25) 



where 7 = m±/m\i . 

Before proceeding to calculation of exchange interac- 
tion between elliptic valleys we first remind several details 
about single-electron ground state that will be important 
for the calculation. Ground state energy spans the inter- 
val changing from —2 to — | when 7 ranges from to 1 
(see Fig. QJ. 7 = 1 corresponds to the spherical sym- 
metric case with the ground state energy Et = — | of a 
hydrogen-like center. For 7 < 1 ground state energy was 
calculated ini^ by variational two-parameter fit with the 
trial function 



I 2 \-l/2 -\/p 2 / a i+ z2 /a» 

(7ra]_aii) ' e V ^ « 



(26) 
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Figure 1: Binding energy of electron in an elliptical valley with 
the transverse to longitudinal mass ratio 7. Solid curve — re- 
sult of numerical computations of the present paper. Dashed 
curve — result of variational fit with the trial function 12611 . 
Solid straight line — approximate adiabatic solution 12711 near 
7 = 0. 



(dashed curve in Fig. QJ. 

Our results of numerical calculations by the method 
of differential sweeping in an infinite interval with 
singularitjii^ are presented in Fig. Q]as solid curve. The 
wave function was sought in the form of the expansion 
in L = 40 first Legendre polynomials in polar angle 9 
with respect to the axis d u . The cut-off at L — 40 of the 
indefinite basis cannot affect noticeably the accuracy of 
calculations until the ratio of characteristic longitudinal 
and transverse length scales 7 1 / 2 becomes comparable or 
less than the minimal scale ~ L^ 1 of the truncated ba- 
sis. For that reason we have not plotted the numerical 
curve ^11(7) down to 7 — * 0, where larger and larger ba- 
sis would be required. However, an approximate solution 
may be found in an adiabatic approach^: 



£11 



-2-4^^27", 



(27) 



where a' Q w —1.02 is the first root of the derivative Ai' 
of the Airy function. This expression is shown in Fig. 
as solid straight line. 

At large distances from impurity the wave function of 
a localized state with E\\ = —a 2 /2 asymptotically falls 
off as 



tp « A(0)r 



■qj a— 1 -or 



, /4 

-R/2 



-®- 
11 



JL 

R/2 



Figure 2: Frame (£, for calculation of exchange between 
two elliptic valleys with the same d. The £-axis connects the 
two centers. 



where r 



D 



-1/2 



Tj and D l} 



-1/2 



(5ij - dkdj) + 



7 1 / 2 didj. We introduced 



7 1 cos 2 1 



(29) 



where 9 is the angle between C and d. 

In real semiconductors several degenerate elliptic val- 
leys with different d„ (d u along the rotation axis of the el- 
lipsoid) contribute, generally speaking, to the wave func- 
tion and hence to the exchange energy. Before addressing 
the complications due to the valley degeneracy, consider 
first two donors in a hypothetic tetragonal semiconduc- 
tor with only one elliptic valley allowed by symmetry. 
We assume again the two-electron wave function to have 
the form 111411 . where, however, the two single-electron 
functions ip VA (jiA), v^^b) should now be found from 
solution of the Schrodinger equation for an elliptic val- 
ley. Proceeding as above, from a bilinear kinetic energy 
Dijdidj we obtain instead of Eq. (II of) 



R/2 - Ci R/2 + C2 



1 1 
1 

R 



(30) 



ri2 



At large distances r w C,C along the ^-direction r ~ 
q{9)\(/\, where q{9)is from l|29|l . and then from 12811 obtain 



Di 3 dje 



1(0) 



(31) 



(28) 



In the coordinate system chosen the angles 9 are equal 
for the two impurity centers (see Fig. [2j. Hence 1(30}) 
reduces to 1)1511 with a — > a/q(9). This substitution also 
gives the new expression for the matrix elements J to 
replace ltTT|l . Also the two-electron wave functions l(Hjl 
in expression ltTT|l now contain single-electron functions 
corresponding to elliptic valleys. 

To calculate J we now follow the scheme presented in 
Section |n] for atoms. To expand the exponents in the 
single-electron wave functions we supplement the inter- 
impurity axis C with perpendicular axes £ and r) so that 
d = C cos 9 + £ sin 9 and expand r in the small ratios £/£ 
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and rj/( up to second order 



«ICI 



1 



sin 9 cos 9 £ 
? ? 



e 



2<Z 2 C 



2/*2 



(32) 

Applying this expansion to the sum r%A + r ib of the ex- 
ponents we substitute ( with (± = ( ± R/2. Since for 
the electrons located between the ions —R/2 < C < -R/2: 
|C±| = ±C±> terms linear in £ in (l32f) cancel each other. 

Hence like in atomic case we are left with the Gaus- 
sian integrals over transverse coordinates p 1 = (£i,r?i) 
and p 2 that decouple from the integral over £. By the 
substitution (l2fil) the Gaussian integrals reduce to a single 
integral over the polar angle <pi2 of p 12 : 



S(9) 



1 

27 



2tt 



[sin 2 y> + cos 2 ip/'yq 2 



(33) 



In the spherically symmetric case j = q = 1 the integral 
S(0) = 1. 

Calculating matrix element ltTT|l gives eventually 



J = -(4TT) 2 ^A 4 (6)S(9)I(a,q)R 7q/2a - 1 e- 



2aqR 



(34) 



where 



I(a,q) 



T(q/2a) q 3 +i/ 2a 



t 3 / 2a (2-t)^ 2a e- qt/a dt. 



22+q/a a 2+q/2a 

(35) 

and the integration variable in 1(35)1 t = 1 — 2(/R. For 
two spherically symmetric (7=1) atoms with angular 
momentum projection m = the asymptotic coefficient 
A{6) = A/y/In (see (PH) and Eq. (gH reduces to {HI 
as it should. 

In addition to finding the binding energy we have also 
carried out numerical calculation of single-electron wave 
functions to obtain asymptotics 12811 and exchange inte- 
grals (134*)) using the differential sweeping method^. The 
preexponents of the exchange integrals, J(R)e 2aqR , are 
presented in Fig. O^for 7 = 0.1, 0.4, 0.7, and 1, and for 
9 = 0, 7r/6, 7r/4, 7r/3, and ir/2 (from top to bottom) 
for each 7. In the interval of inter-center separations 
1 < R < 10 accessible for numerical evaluation the de- 
pendence of preexponents on R is approximately linear 
in double logarithmic scale with the slope agreeing well 
with 7q/2a— 1 which follows from H34H . 

Returning now to realistic semiconductors, recall that 
degenerate valleys have different anisotropic dependence 
for the decay of the wave function of localized electron. 
For an arbitrary orientation of the axis connecting two 
donor centers the main contribution to exchange inter- 
action comes from the wave function components corre- 
sponding to the pair of valleys with the slowest decay, 
contributions of the other pairs being exponentially sup- 
pressed. In other words, among all the elements 

(uv\J\u'v') ^e- Q ^") + ^ e ") + ^ e "') + «^')] fl / 2 , (36) 

we can, with the accuracy up to exponentially small 
terms, keep only the elements between such states u that 




Figure 3: Preexponents of exchange integrals between ground 
states of electrons localized in elliptic valleys vs. donor center 
separation for four values of 7 = 0.1, 0.4, 0.7, and 1. For 
each value of 7 five curves are plotted for different angles 
between the preferred direction of the elliptic valley and the 
inter-center axis: 9 = 0, vr/6, 7r/4, tt/3, n/2 (from top to 
bottom). In the spherical symmetric case 7 = 1 exchange 
does not depend on 9. 



the exponent in l-'Uill is minimal: (uu\ J\uu) ~ g-2a<j(e u )-R_ 
Therefore, for any arrangement of impurity centers ex- 
cept for high symmetry directions the electron exchange 
is predominantly caused by the overlap of the wave func- 
tion components from valleys with cIa = (1b, which 
means that the term splitting is correctly described by 

Eq. eg. 

The relative position of the two donors being changed, 
crossover between the pairs of valleys contributing the 
most into the exchange matrix element ll-'ifit occurs 2 ^. 

These considerations are not changed if one takes into 
account the lifting of the inter- valley degeneracy of a sin- 
gle center by the short-range part of the impurity poten- 
tial ("central-cell correction"). The ground state in this 
case is the non-degenerate symmetric combination of the 
functions QJ l0 



r 



1 



v U — 1 



(37) 
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where ip u { r ) oc e lk " r are the Bloch functions at the mini- 
mum of the conduction band located at k = k u , and F is 
the number of the participating valleys (6 in Si, 4 in Ge, 
and 3 in III-V compounds). 

Considering the Bloch functions oscillate rapidly on 
atomic scale much smaller than the effective Bohr radius, 
exchange interactions between the states (|37l may be 
seen to be 3 - 



1 



J 9 « ^ E (uv\J\uv)e^ 



k„)R 



(38) 



u,v—l 



If we average over the rapidly oscillating exponential 
terms, we get 



a single-electron state. The permutation operator P in 
© now interchanges both spin and coordinate variables 
of the electrons. 

Correspondingly, the provision for the total wave func- 
tion ■0 © to be antisymmetric in electron permutation: 

(p) 

ip = -Pip, leads to c„ A „ B = -ci, AVB . 

For simplicity we will restrict ourselves to donors in 
II- VI and III-V semiconductors where besides being a 
Kramers doublet single-electron state on a donor is not 
otherwise degenerate. Then v will stand for the Kramers 
index. 

In bulk zinc-blende semiconductorsiii^ 



Mk) = 



zk x (k 



(41) 



I L 

u=l 



(39) 



which shows that the exchange between the symmetric 
ground states of multi-valley donors for any arrange- 
ment except along high symmetry directions is described 
by Eq. II34J1 with the asymptotic parameter A(&) — > 
A(9)/VT. 

Concluding this Section, we discuss briefly the pre- 
vious study of exchange interaction between donors 3 . 
In the cited paper, the anisotropy of the exponential 
dependence of exchange interaction for elliptic valleys 
was completely ignored by substituting the expression 
for isotropic atomic exchange. Whereas, e.g., for Si 
7 = m±/m\\ ~ 0.21, and hence q(6) Q29JI in the ex- 
ponent of the exchange interaction oc e 2aq ^ R changes 
from 7- 1 / 2 « 2.2 for 6 = to for 9 = tt/2. Besides, 
in Ref»i the Heitler-London approximation was used, the 
approach proved incorrect in calculating the preexponen- 
tial factors at large distances. On the other hand, we ob- 
tained exact asymptotic expressions for the dependence 
of the exchange interaction on the valley parameter 7 and 
the mutual orientation of donors 8. 



IV. ASYMMETRIC EXCHANGE 
INTERACTION (DONORS) 

In semiconductors without center of inversion symme- 
try, such as III-V or II- VI compounds, spin-orbital inter- 
action for donors acquires the for m 1 ^^ 



Hso = h(-id)s. 



(40) 



Wave function of a single electron cannot be represented 
as a product of spin and coordinate components any 
more. An approach with the expansion 10 is nevertheless 
applicable after suitable adjustments^. 

In the presence of spin-orbital Hamiltonian 14011 spin 
projection no longer commutes with the total Hamil- 
tonian, but the two-fold (Kramers) degeneracy of the 
ground state remains. Index numbering this degeneracy 
may be included into the full set v of indices describing 



and h y (k) and h z (k) are obtained by permutation in l|41H . 
Here E g is the band gap in atomic units and a so is a 
phenomenological parameter. For GaAs E g » 1.43 eV 
(in atomic units this yields E g w 113 under the radical 
sign in H4H ) and a so « 0.07. 

For small spin-orbital interaction the spinor single- 
electron wave functions are 



<^(r) ~ <fo{r){e 



tot h(f)ro r /2\/3f 



(42) 

where <po(r) — e~ r j \px is the hydrogenic ground state 
wavefunction and f = r/r. The phase factor turns into 
unity when r — > 0. Thus the Kramers index takes on 
the meaning of the spin projection of an electron on the 
center v =f , \ . 

Two-electron wave functions have now to be sought as 
solutions of the Hamiltonian 

H = -±dl-±d 2 2 +H S0l + H S0 2-l/r 1A -l/r2B 
- 1/raA - 1/riB + IA12 + 1/R. (43) 

in the form II17|I 
<U(*i,*2) = * aa '' ^>i,r 2 )<(r 1A )^>2B), (44) 

with x aa converting to 5 aa <5 /3/3 when either Tia — ► 
or r 2 e — > 0. The principal terms in spin-orbital interac- 
tion are included into definition (@2j, and the equation 
for x aa turns out to coincide with 11511 , Therefore 

x aa\/3/3' = ^aa'^' with % from ^ 

Similarly to iflUjl we get 



(E - E )c v 



V J , , c {P) 



(45) 



where 



= J{R){e-^' 2 ) VA ^{e^l 2 )^ A . (46) 
The exchange integral J(R) = -0.818R 5/2 e~ 2R £& and 
7 = HC)R- (47) 



8 



Absolute value |"y|/2 plays the role of angle of rotation of 
electron spin around the axis 7 when it tunnels between 
each of the centers A and B and the median plane. Using 



(48) 

we finally arrive at the equation [E — Eq)c = 7i ex c for 
the coefficients c VAVB , with the HamiltonianiS 



(49) 



<ti R (7)6-2 + 1 



PS -ij [0-10-2 + 1 - 7 (0-1 x 6- 2 )] , 

where the Pauli matrices &\ and <x 2 act in the space of 

< — > 

the Kramers indices, and R (7) is the three-dimensional 
rotation matrix on the angle 7 around the axis 7. 

In the absence of spin-orbital interaction 14011 the angle 
I7I = and we return to the exchange Hamiltonian <j£| 
as expected. 

In this Section we have considered the expansion of 
spin-orbital interaction 14(H in bulk III-V semiconduc- 
tors, where ijlTijl begins with the third power in momen- 
tum. Analogously, one can study spin-orbital interaction 
linear in momentum which appears in reduced symme- 
try or dimensionality (in films or at the surface). All the 
derivation including the expression (HTjl for angle 7 re- 
mains intact. The case of linear spin-orbital interaction 
has also been studied differently in Refi^fi. 



V. DEGENERATE BANDS (ACCEPTORS) 

As we will see later, calculation of the exchange inter- 
action between two shallow acceptor centers has much in 
common with the problem of two hydrogen atoms each in 
a state with non-zero momentum discussed in Section ITT1 
and in the Appendix. However, the analogy is not com- 
plete, the calculations become laborious and we have not 
performed them fully. In what follows we only determine 
the sign of the exchange interaction. More precisely, it 
is shown that the magnetic moment in the ground state 
of a pair of acceptors equals zero. 

We begin this Section by first reviewing the properties 
of a single acceptor center. 

Without taking spin-orbital interaction into account 
the valent band of typical semiconductors (Ge, Si, III-V) 
would be six- fold degenerate at the point k = 0. Spin- 
orbital interaction splits the valent band into a two- and 
four-fold degenerate bands separated by a gap. If the 
spin-orbital splitting is large with respect to the binding 
energy of an acceptor, the two-fold degenerate band can 
be neglected. For k ^ the four-fold degeneracy of the 
upper-lying band is lifted and two two-fold bands of light 
and heavy holes are formed. Their energy spectrum e(k) 
quadratic in k is known to depend on the direction of 
the quasimomentum k through three invariants of the 



Oh point symmetry group so that surfaces e(k) = const 
are slightly warped spheres. 

As is known, the spectrum in the vicinity of k = 
is determined by three Luttinger parameters 71, 72 and 
73. In practice, two of them are close 72 ps 73 = 7 and 
the dispersion relation is simply e = k 2 /2mi l h, where 
fni,h — m o/("fi ± 7) are the effective masses of the light 
and heavy holes. In this so called spherical approxima- 
tion the Schrodinger equation for the four-component en- 
velope hole wave function tp a is given by Q with the 
effective mass tensor 



Dij = — [71* 
m 



y -7(Wi}-|J%) 



(50) 



where J is the 4x4 matrix of the pseudospin operator 
J = § and {JiJj} — ^(JiJj + JjJi)- 

It is convenient to rewrite 15011 using atomic units <2b = 
At/m^e 2 and 2Ry = l/rrihC%: 



Di 



Here 



_ 1_ _ mh - mi 
^ 7i m h +mi' 



(51) 



(52) 



In the limit of equal masses m/, and to/ the four 
components of ip a in (2J decouple, each satisfying the 
Schrodinger equation for a hydrogen atom with the 
ground state energy Eh = —1/2 for orbital momentum 
L = 0. In the general case mu > mi acceptor states 
may be classifie d 21 i 22 ' 2 ^ by the total angular momentum 
F = J + L and its projection m F : 

\Fm F ) = ^ X! (mjm L \Fm F )\Jm J )\Lm L } 1 (53) 

L mj,m L 

where (mjm^FmF) are the Clebsh-Gordan coefficients 
and mj + m L = m F , \Lm L ) = R L (r)Y£ lL (r), Y F riL (r) 
are the spherical harmonics and \ Jmj) is a column in 
the pseudospin index with one non-zero component at 
to j. The summation over L for a given F runs over all 
such values that \L — J\ < F < L + J. 

Substituting ffify into © with <OJ results in the 
Sturm-Liouville problem on the energy Eh and the un- 
known radial functions Rl{t). The ground state turns 
out to be four-fold degenerate with F = | (It is, in fact, 
the only state with possible L = 0, hence it becomes the 
correct ground state in the limit to^ = to;). In this state 
only the radial functions Rl{t) with L — and 2 turn 
out to be coupled: 



-M (d 2 r - \d r 



-li (d 2 r + &d r + 4r) 

dl + la, - « 



+2(1-/1) + 



Ro 
R2 
Ro 
Ri 



(54) 
= 



Numerical solution 2 ^ of this equation shows that E h 
changes from — | at /1 = (to/j = mi) to ps — | in the 
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limit /i — > 1 (rrih 3> mi). Thus the binding energy is 
determined by the scale of the heavy mass. At large r 
|H|| yields 



(l-/i)(3 r 2 + 2^)(i? 
(d? + 2E h mi/m h )(Ro 



J2a) 
i? 2 ) 



(55) 
(56) 



which shows that for to/j ^> mi the sum of the two func 
tions falls off on the scale 1 



2\Eh\ corresponding to the 
heavy holes, i.e. much faster than their difference, falling 
off on the scale determined by the light holes. Then at 



large 



■24 



Ro 



— i?2 oc e 



with 



cti = y / 2\E h \m l /m h . 



(57) 



The four degenerate wave functions of the ground state 
l|53l) in the asymptotic region close to the axis £ connect- 
ing the impurities (chosen as the quantization axis for 
the orbital momentum L) are 




(58) 



Pf(r) 
</>|(r) 
V-i(r) 

</?-|( r ) 



Here the four unit columns \ Jmj) = | ± 3/2), | ± 1/2) 
form the basis in the pseudo-spin 3/2 space, and we have 
left only terms up to first order in the small polar angle 
6 (cf. Eq. CHI). 

Proceeding now to the calculation of exchange inter- 
action, we construct asymptotically correct basis of the 
direct product of the ground states on two centers in the 
form 



a/3 



1 (r 1 ,r 2 )= X a7 ' /35 (ri,r 2 )^. 



(riA)v m B(r 2B )- 



(59) 

The total two-electron wave function ip is again sought 
in the form (HjJI, where the condition c m A TO B = — c I B 

is required for the wave function tp to be antisymmetric 
in particle permutation. Equation on \ is obtained by 
substitution of this form into the two-hole Hamiltonian: 



1 



1 



D 



1 



R/2 - Ci R/2 + C2 R 



T\2 



a-y,/35 



= m 



where we used the asymptotics dup f 
-QWpQ sgnC at large r = CICI (cf- ©30- 

We neglect the derivatives of x over transverse coordi- 
nates as they are, as before, of higher order in 1/R (the 
validity of this approximation is verified by differentiat- 
ing the resulting expression for \). 

Moreover, in the spherical approximation we consider, 
the hole Hamiltonian is invariant under transformations 



from the spherical point symmetry group Kh and not 
only from Oh- Therefore, in this approximation, axis C 
connecting the acceptor centers may be chosen as the 
quantization axis of the moments J, L and F, and the 
problem acquires some resemblance with the hydrogen 
problem studied in Section I. 

In particular, the equations ll6f))l on x decouple from 
each other: 



ai(D«fd ClX a '^ S - D^'d C2 x a ^' S ) + 



1 



1 



1 



R/2 - Ci R/2 + ( 2 R 



1 

r\2 



X 



efy,f38 



$51) 



because matrices D"" and are diagonal: 



D zz = 



( 1 

mh/mi 

mh/mi 

\ 1 



as can be easily seen using the explicit form of the pseudo- 
spin matrix J z . Together with the boundary condition 
% "7,/3<5 s^S? 6 when either ri -> -R/2 or r 2 -> R/2, 
this equation yields solution for x" 7 '* 3 * diagonal in each 
pair of band indices, i.e. only the elements with a = 7, 
(3 = 5 are non-zero. 

Because of such diagonal structure, the one-hole func- 
tions in the elements of the exchange integral matrix 



J, 



1 A l y~ l - l ES fy-y, A 



a/3 " / &-6> 



(62) 



x <a (i-ia)^b (r 2B )^A, (i"2a)^b , (r 1B )dCd i p 12 . 

At integrating wave functions II58J1 over transverse co- 
ordinates in l(62l) each 9 will produce an additional small- 
ness RT 1 ! 2 as we have seen in Section I. In princi- 
pal order we may hence leave only functions ip ± i{v) « 

Rq (r) I ± / (each of them has also only one non-zero 
pseudospin component). 

The 16x16 diagonal matrices D®" S 1313 before and 
fiaa £)/3f3 j-, e f ore m do no t generally coincide for 
a / j3, and hence equations for all the elements of X 
cannot be simply solved in terms of C12 as was done in 
atomic case. Nevertheless, to calculate II62|) in main order 
we only need ip±i as we explained above. Equations 16 111 
on x for the components ip± 1 do reduce to l|15ll with a — > 
aiirih/mi. We thus arrive to a complete analogy with the 
problem of exchange between two hydrogenic centers in 
the ground state, i.e. degenerate only over projection 
of spin h. In this analogy localized wave functions are 
Ro(r)/y/n and the role of the projection of "spin" | is 
played by to_f spanning the truncated subspace ±|. 

We may, therefore, at once give the asymptotically cor- 
rect terms: the ground state will be the "singlet" 



V2 



(63) 



10 



and the excited state will be the "triplet" 



V2 V 



I _I\ 

2 ' 2/ 



(64) 



where -0 = xRo( r iA)Ro(r2B)/^ and x is given by ifTTjl 
with a = aimh/mi. 

Energies of the "singlet" and "triplet" states are Eq ± 
J(i?) respectively, where for the exchange splitting we 
obtain 



the region of localization r ~ ab, and Eh ~ — \- High 
magnetic fields ob S> an localize an electron in the trans- 
verse directions much stronger than Coulomb potential 
does, and the binding energy grows logarithmically with 
magnetic field Eh — —\ In (ae/a//) 2 . 

Exponential tails of localized electron wave function in 
magnetic field have the following form (see^ and refer- 
ences therein). For a weak field (an > a^) there are two 
asymptotic regions with different dependencies of ip(r). 



Rather close to the z 
functions fall off as 



axis, 1 <C r <C (ajj/os) , wave 



mi 
m h 



1 

2^ 



2 i? 3 



-2aR 



if oc exp —r 



1 + 



24 



(67) 



f 1 -P 2 [f (1 + 0] -P 2 (1 t)] n Lzl^ (f ,f.s whereas at large r > (a^/as) 2 as 

i ri + ti*-»fi-t)i-» ' ( } 



where P(r) is defined as P(r) = R(r)e 2aR with R(r) 
from numerical solution. For two atoms, when mi — mh, 
on = a and the preexponent P{f) = Ar 1 / Q_1 /2, Eq. (I f i o f) 
reduces to IpU- 

Asymptotics Ro(r) w -R2M — R(r) used to obtain 
J(i?) unfortunately sets in too far (i? > 10) and we have 
not evaluated Q6FiJI numerically. 

Similar to the hydrogen molecule problem the two 
states found, 16-tll and Il(i4t . with energies split by the 
amount of 2J(R) l|65|) represent only the ground and 
the most excited states. Energies of all the other two- 
center states are defined by smaller exchange integrals 
and will have exchange splitting 0(R^ 1 J(R)), i.e. at 
large R these terms will lie between the first two. 

As for the magnetic moment of a single acceptor center, 
it is proportional to F with the g-factos^ g « 0.76. We 
see immediately that the ground state of two acceptor 
centers is non-magnetic as is the case with hydrogenic 
atoms and donors. 



if oc exp - 



1*1 + 



2 / \ 2 
P ( a B 



4 \a H 



(68) 



In strong magnetic fields (ajj Cob) the first asymptotic 
region is absent and exponential tail of the wave function 
is given by 



ip ex exp - 



Iln(^ 



EL 

4 



flB 



(69) 



The change in the "longitudinal" asymptotics (field along 
the direction connecting two centers) compared to l|68jl 
is due to the change of the binding energy (and hence, of 
the localization length) with magnetic field. 

From these formulas we observe that application of 
magnetic field parallel to the inter-center axis does not 
change the form of the exponential dependence of ex- 
change integrals. At the same time, field perpendicu- 
lar to the axis changes exponential asymptotics either to 
slightly corrected atomic one 



VI. ROLE OF MAGNETIC FIELD 



,-2R(l+R 2 (a B /a H r/6) 



(70) 



To complete the above discussion, we briefly review 
the impact of external magnetic fields. Strong magnetic 
field squeezes electron wave functions in cross direction 
reducing their overlap. Therefore, major changes come 
already in exponential factors, and one may use well- 
known results for the asymptotic behavior of hydrogenic 
wave functions in magnetic fields. 

Uniform magnetic field H = Hz, with the vector po- 
tential A(r) = x r results in an additional magnetic 
parabolic potential in the Schrodinger equation for the 
ground state 



2 r 8 \a H 



(66) 



where a 2 H = hc/\e\H is the magnetic length. In weak 
magnetic fields ob <C clh magnetic potential is weak in 



for 1 <C R <C {aH/a 2 B ) in weak fields, or to 



-R 2 (a B /a H ) 2 



(71) 



for R » (an/ag) in weak fields and for all R S> 1 in 
strong fields. 

We conclude that magnetic field drastically changes 
already the exponential dependencies, and we thus will 
not study the behavior of the pre-exponential factors. 



VII. SUMMARY 

We have calculated exact asymptotic form of exchange 
interaction between shallow impurity centers in semicon- 
ductors at large separations generalizing the methods of 
the theory of hydrogenic molecules. 
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We found that the ground state of a pair of impurity 
centers is always non-magnetic. In other words, the ex- 
change interactions between centers have the "antiferro- 
magnetic" sign. 

The essence of the method that we use is the construc- 
tion of a secular equation with the matrix comprising ex- 
change integrals between all the degenerate single-center 
states. The exchange integrals are expressed as the prob- 
ability current through the hyperplane in the coordinate 
space of two electrons that separates two regions of their 
localization on the centers. The exchange integrals are 
calculated with wave functions at the hyperplane cor- 
rected for the Coulomb interaction between electrons and 
alien ions in principal order in inverse inter-center sepa- 
ration R . 

Exchange integrals between different degenerate single- 
center states may have different order in R~ . The main 
exchange integrals define the energy of the lowest and 
the highest split terms with the remaining terms lying in 
between. To calculate the latter the exchange integrals 
in the next approximation would be necessary. However, 
the procedure simplifies for atoms if one uses symmetry 
considerations to restrict the basis of degenerate states on 
each center to only those that are allowed by symmetry 
to convert to a given term. 

In semiconductors we have considered degeneracy of 
single-center states both in valley index (donors) and in 
the band spectrum near k = (acceptors) . For electrons 
localized on donors calculation of exchange interaction 
reduces to evaluation of exchange integrals between the 
components of the wave functions from different pairs 
of valleys. Depending on the relative position of the two 
donors in the cubic lattice, crossover from one pair of val- 
leys to another takes place at high-symmetry directions. 
The dependence of exchange integrals on valley param- 
eters and on the orientation of the valleys with respect 
to the axis connecting the two centers was calculated nu- 
merically. 

We have also evaluated the asymptotic form of the 
asymmetric (Dzyaloshinskii-Moriya) exchange interac- 
tion for donors. 

The problem of exchange interaction between two ac- 
ceptors simplifies considerably in the spherical approx- 
imation when it acquires resemblance with hydrogenic 
molecules. We have developed the complete procedure 
for calculation of exchange in spherical approximation, 
but have restricted ourselves to the sign of the interac- 
tion only. As in the case of a hydrogen molecule, the 
ground state of a pair of acceptors turned out to be non- 
magnetic. For a detailed calculation of the exchange in- 
teraction numerical solution of the wave equation of an 
acceptor center at very large distances would be required, 
which was not possible to do in the numerical scheme we 
used. 

The influence of external magnetic field on the strength 
of exchange interaction was also briefly discussed. Ex- 
pectedly, magnetic field affects already the exponential 
dependence due to stronger localization of electron on a 



center in the presence of magnetic field. 
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Appendix A: HIERARCHIC STRUCTURE OF 
MOLECULAR TERMS 

In this Appendix we present asymptotic expressions 
for exchange integrals between higher orbital momentum 
I atomic states. Preexponents of these integrals contain 
R in different degrees dependent on the projections m of 
I for the two atoms. We show how exchange integrals 
(and hence electronic terms) of higher than leading order 
in R^ 1 may be found by restricting the basis of single- 
atom wave functions to only those allowed by symmetry 
to participate in a given molecular term. 

The matrix elements (1 1 H for arbitrary I were consid- 
ered in Refiii in the most general case of different bind- 
ing energies of the two atoms. However, these general 
formulas were cumbersome and were not brought to ex- 
plicit answe»2I. Below we cite our results for the case of 
equal binding energies: 

x R^- l -^e- 2aR , (Al) 

where 

r(l/2a) , . , 

I a ± = — r— ™+ K{2a,mi,m 2 , Am) 

X / (l_y)£±— e ^ dy ( A2 ) 

Jo 

and 

m+ = \rriA\ + \m' A \ + |m B | + |m' B | 
m- = |toa| + \m' A \ - (|tob| + \m/ B \) , 

mi = \m A \ + \m B \, 

m 2 = \m' A \ + |m B |, 

Am = \m A - m' B \. (A3) 
The four-fold integral 

K (2a, mi, m 2 , Am) = ^1/2^ (1/ 2a) J ^VV™ 2 
xe -p?-P2 cos [Am(ipi - ip 2 )] d 2 p 1 d 2 p 2 (A4) 




Figure 4: Constants Ji(ct) and J2(«) from Eq. llA7(l . 



may be calculated analytically. First three integrals are 



K(2a, 0,0,0) = 1, 
#(2a,2,0,0) = l + l/4a, 
A"(2a, 1,1,1) = -l/4a. 



(A5) 



(See alsoii). 

We observe from Eq. IjAlll that the matrix elements for 
transitions between states with different m have different 
degrees of R in the preexponents. Since in all the calcu- 
lations we have omitted all but the leading terms in R , 
secular equation ltl3|l is valid only up to the accuracy of 
the matrix element with the highest preexponent. 

In the complete basis of all the degenerate atomic 
states the principal matrix element in is the one 
with minimal m + = 0, i.e. between the states with 
771 a = 777b = m' A = m' B = l|21|) . It corresponds to 
one of the molecular X states. With the accuracy up to 
this element, matrix J mA m B m' m' is trivial with only one 
finite element Jq{R) at ttia = tt7b = Tn' A = m' B = 0. So, 
at large R molecular terms are arranged as follows: the 
ground singlet E state with the energy E Q — \J (R)\, the 
triplet excited £ term E + \J (R)\, and the other terms 
with energy is E Q + 0(R^ 1 J a (R)). 

One can resolve the energies of those terms by restrict- 
ing the basis of degenerate atomic functions to only those 
that by symmetry can convert to the molecular term with 

given A = |TOA + 777b | = |777 A + 777^ |. 

Because of the angular momentum projection conser- 
vation 777A + 777b = 777 A + 777g the next minimal 777+ = 2 
will be for the matrix element between the states with 
one of 777a, itt-b — 1, the other zero. Such atomic states 
(and only them) convert to molecular II states, which 
thus will have splitting Ji 7 / 2a ~ 2 e ~ 2a R^ Arranging the 



four states in the following order |777at77b) = |0, 1), |1,0), 

E-En 



\Jl\+Jl 



\Jl\-J-2 



-W + J2 



singlet |01) - [10> + P\01) - P|10) 

triplet |0l) + |10) - P|01) - P|10) 
singlet |0l) + |10) + P|01) + P|10) 



\Ji\-J 2 -- triplet |0l) - |10) - P|01) +P|10) 



Figure 5: Sketch of the energy level splitting of p-states 
of hydrogenic atoms under exchange interaction in diatomic 
molecule (in units of ±|G| 4 (2Z + l) 2 l(l + l)R 7 / 2a ~ 2 e~ 2aR ). 



|0, — 1) , I - 1,0), we find that the matrix J mKm nm' A m' B 
J = -\Ci\ i {2l + lfl{l + l)R 1 l 2a - 2 e~ 2aR 



is 



/ J 2 Ji 

Ji J 2 

J 2 J 1 

V Ji J 2 



(A6) 



where 
Ji 
Ji 



-1(2 + J- 

r(l/2a) 



7(a), 



2 3 +ia 4 +^ 
1 



"(1 - y)^-\l + y)^-\l + y 2 )diA7) 



These constants are plotted as functions of a in Fig. 21 

Matrix <l AB|) has two doubly-degenerate eigenvalues 
J 2 ± Ji = J 2 T l^i I- Hence eight-fold degenerate two- 
electron state will split into four doubly-degenerate co- 
ordinate states located in the order shown in Fig. El 
The correct combination of two-electron wave functions 
converting into molecular terms are presented in Fig. 
for the total projection of the orbital momentum L = 1. 
The second degenerate state with L = — 1 for each term 
is given by substitution of 1 with —1 in the two-electron 
wave functions. 

Note that unlike the energy level splitting of s-terms, 
the lowest energy term of the two atoms each in a p- 
state turns out to be spin-triplet while the most excited 
— spin-singlet. 
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